基于 高 空 化 性 能 的 对 转 泵 喷 推 进 器 设计 “ 


胡 雷 俊 "， 曹 琳 琳 * ， 车 邦 祥 ， 吴 大 转 


(1 浙江 大 学 能 源 工程 学 院 化 工 机 械 研究 所 ， 浙 江 杭州 310027; 
2 浙江 大 学 流体 动力 与 机 电 系 统 国 家 重点 实验 室 ， 浙 江 杭州 310027) 


摘 要 : 推进 器 将 动力 设备 发 出 的 机 械 能 转化 成 推力 以 驱动 航行 体 前 进 ， 是 决定 航行 体 性 能 的 关键 动力 装置 。 随 航速 提升 ， 

进 器 易 遭 遇 空 化 、 品 声 与 推进 效率 下 降 等 问题 。 结 合 高 速 推 进 的 需求 以 及 对 转 泵 的 性 能 优势 ， 本 文 提 出 对 转 硝 喷 推 进 器 结构 ， 
通过 前 后 转子 串联 、 反 向 旋转 的 形式 ， 以 更 高 密度 的 能 量 转化 模型 蔡 代 传统 模式 ， 在 同等 设计 参数 下 实现 转子 转速 、 叶 片 载荷 
均 相 对 较 低 ， 从 而 提升 航行 体 推 进 器 的 抗 空 化 性 能 及 相关 声学 性 能 。 文 中 基于 某 航行 体 性 能 参数 ， 通 过 相似 换算 设计 了 实验 模 
型 泵 喷 及 其 实验 系统 ， 并 基于 数值 模拟 对 推进 器 性 能 进行 了 研究 。 结 果 表 明 ， 对 转 泵 喷 推 进 器 叶片 的 欧 拉 能 量 分 布 及 载荷 模型 
明显 优 于 单 转子 推进 器 叶片 ， 抗 空 化 性 能 也 有 明显 提高 。 
关键 词 ， 泵 喷 推 进 器 ， 水 下 航行 体 ， 对 转 泵 ; 抗 空 化 设计 ;叶片 载荷 分 布 
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Design of Contra-Rotating Water jet Pump Based on 


High Cavitation Performances 


HU Lei-Jun', CAO Lin-Lin'', CHE Bang-Xiang!, WU Da-Zhuan'? 
(1 Institute of Process Equipment, College of Energy Engineering, Zhejiang University, Hangzhou 310027, China; 
2 State Key Laboratory of Fluid Power & Mechatronic Systems, Zhejiang University, Hangzhou 310027, China) 


Abstract: The propulsor converts the power from engine into thrust to push autonomous underwater vehicle forward. 
With the increase of advance speed, the propulsor is possible to suffer with cavitation, deterioration of acoustics 
performances and propulsion efficiency. In this paper, the contra-rotating waterjet pump with two rotors rotating 
oppositely is introduced as a kind of propulsion device with high density of energy. Under same design specification, 
the rotation speed, blade loading as well as cavitation performances are supposed to be superior to the traditional 
single rotor type. Based on the parameters of certain vehicle, the experimental model and system are designed in this 
paper. And the loading status and Euler energy distribution of impeller blade in contra-rotating type are distinctively 
more reasonable than the single rotor type, and the cavitation performances are obviously improved. 

Key words: Waterjet pump; Autonomous underwater vehicle; Contra-rotating pump; High cavitation performances 
design; Blade pressure distribution 


0 引言 需 推力 及 工作 下 洪 深度 等 因素 直接 决定 着 推进 系统 
的 形式 。 现 在 航行 体 主流 的 推进 器 形式 为 单 转子 螺 

推进 器 将 动力 设备 发 出 的 机 械 能 转化 成 推力 旋 桨 、 对 转 螺旋 桨 、 泵 喷 推 进 器 及 导管 螺旋 桨 。 从 

以 驱动 水 下 航行 体 前 进 ， 因 此 航行 体 航 行 速度 、 所 技术 角度 来 看 ， 无 论 是 螺旋 桨 还 是 泵 喷 推 进 器 ， 随 
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着 航行 体 航 速 与 推进 功率 的 提升 ， 其 推进 系统 通常 1.1 对 转 泵 喷 推 进 器 参数 


遭遇 一 系列 高 转速 带 来 的 问题 ， 其 中 最 核心 的 即 空 本 文 基于 茶水 下 航行 体 性 能 参数 ， 通 过 相似 换 
化 、 噪 声 和 推进 效率 1 I, 算 设 计 了 第 一 代 实 验 模型 泵 喷 推 进 器 。 在 最 初 设计 


传统 形式 的 推进 器 提升 推进 力 多 依靠 提升 转 验证 阶段 ， 第 一 代 实 验 模型 泵 喷 设 计 为 轴 流 转子 形 
速 来 实现 。 随 推进 器 转速 增高 ， 流 体 相对 速度 大 幅 式 ， 由 前 转子 、 后 转子 及 喷 管 三 部 分 构成 。 计 算 模 
提升 ， 高 转速 推进 器 的 空 化 现象 很 难 避 免 。 就 叶片 型 如 图 1 所 示 。 

空 泡 本 身 而 言 ， 其 生长 及 脱落 伴随 有 剧烈 的 宽频 破 图 中 水 流 方向 为 从 右 向 左 ， 从 入 口 向 出 口 看 ， 
裂 噪声 ， 严 重 削 减 推进 器 声学 性 能 四， 空 泡 的 生长 前 转子 沿 逆 时 针 方 向 旋转 ， 后 转子 沿 顺 时 针 方 向 放 
及 脱落 不 断 改 变 着 叶片 表面 压力 分 布 ， 导 致 叶片 非 转 。 计 算 域 入 口 距 离 前 转子 前 缘 约 两 倍 直径 长 度 ， 
定常 载荷 增 大 ， 高 速 旋转 中 各 叶片 空 泡 发 展 程度 又 出 口 距离 导 流 锥 尾部 约 三 倍 直径 长 度 ， 此 设置 已 经 
不 尽 相 同 ， 致 使 推进 器 遭遇 强烈 振动 ; 抗 空 化 性 能 通过 网 格 无 关 性 验证 。 采 用 ANSYS TurboGrid 16.0 
较 弱 的 推进 器 叶片 ， 甚 至 遭遇 空 泡 规 模 陡 增 ， 致 使 对 旋转 体 进行 网 格 划 分 ， 为 方便 生成 结构 化 网 格 ， 
流 道 大 范围 堵塞 ， 做 功能 力 及 推进 效率 急剧 下 降 。 在 导 流 锥 尾部 以 后 增加 了 一 个 极 细 的 虚拟 轮 载 ， 对 
因而 ， 在 保证 同等 推进 能 力 下 转速 较 低 是 推进 器 向 于 水 力 性 能 计算 结果 的 影响 可 忽略 。 

高 速 推进 方向 发 展 时 的 一 个 重要 思路 。 

对 转 泵 是 一 种 由 前 后 两 个 对 旋转 子 串联 而 成 
的 水 泵 ， 其 核心 理念 在 于 前 转子 的 出 流 残 余 能 量 被 
后 转子 充分 利用 ， 转 化 成 压力 水 头 ， 实 现 水 泵 整体 
无 旋 出 流 ， 大 幅度 增加 了 其 能 量 密度 。 后 转子 对 扬 


程 的 贡献 远 远 超出 传统 的 后 置 定子 ， 并 可 通过 转速 — Y — 
的 调节 来 调整 扬程 以 及 扭矩 在 两 个 转子 间 的 分 配 。 (a) 全 流 道 计算 模型 ( 带 泵 壳 ) (b) 全 流 道 计算 模型 (不 带 泵 壳 ) 
通过 充分 利用 前 转子 的 出 流 能 量 ， 相 对 于 同等 设计 (a) Computational model (b) Computational model 
E 数 (流量 ， 扬 程 ， 比 转速 等 ) 下 的 传统 水 和 泵 ， 此 类 of the whole flow passage of the whole flow passage 
水 泵 在 相关 理论 及 实验 研究 中 均 被 证 实 具有 更 低 的 (with shroud) (without shroud) 
转速 、 更 高 的 水 力 效率 ， 更 紧凑 的 体积 结构 、 更 稳 图 1 对 转 泰 喷 推 进 器 计算 模型 
定 的 性 能 曲线 以 及 更 优 的 空 化 性 能 5 o Fig. 1 Computational model of contra-rotating 

对 转 泵 的 前 述 优 势 与 高 速水 下 航行 体 对 推进 waterjet pump 


系统 的 需求 十 分 吻合 ， 虽 然 在 外 观 及 内 部 流动 机 理 
上 上 ， 这 种 对 转 泵 与 安装 导管 的 对 转 螺旋 桨 有 很 多 相 


du Ew Al. 2B pn. 


似 之 处 ， 但 从 推力 产生 机 理 上 来 看 ， 对 转 式 有 泵 喷 依 表 1 对 转 泵 喷 推 进 器 性 能 设计 参数 
然 属 于 水 泵 范畴 ， 利 用 转子 叶片 扬 水 、 流 道 变化 来 Table 1 Performance design parameters of contra-rotating 
实现 系统 高 速水 流 喷射 以 提供 推进 力 。 waterjet pump 

本 文 的 主要 内 容 即 是 提出 将 对 转 水 泵 理念 引 扬程 ”流量 ”转速 比 转速 推力 
入 泵 喷 推 进 系统 ， 形 成 面向 小 型 水 下 航行 体 的 对 转 Him  Q/m's! n/min! — nm, T kN 
硝 喷 推进 器 概念 。 并 基于 某 水 下 航行 体 性 能 参数 ， 前 转子 72 0212 2000 7647 — 
通过 相似 换算 设计 了 实验 模型 泵 喷 及 其 实验 系统 ， 后 转子 72 0212 200 7647 
并 借助 数值 计算 技术 对 其 进行 了 前 期 研究 。 泵 整体 144 0212 s L5 


1 对 转 泵 喷 推 进 器 设计 表 2 对 转 泵 喷 转 子 几何 参数 


Table 2 Geometrical parameters of contra-rotating 


waterjet pump 


HA WRA HMR AA 
D,/mm D,/mm tc /mm Z 
前 转子 200 104 0.5 5 
后 转子 200 104 0.5 6 
1.2 对 转 泵 喷 转 子 设计 方法 
泵 喷 推 进 器 的 高 空 化 性 能 设计 从 两 方面 入 手 : 


一 方面 ， 采 用 对 转 结构 减 小 单 极 转子 的 叶片 负荷 
降低 转子 旋转 的 线 速 度 ， 男 一 方面 ， 在 对 转 结构 
基础 上 ， 进 一 步 对 叶片 豆 型 进行 载荷 优化 设计 。 

对 旋转 子 采用 等 速 差异 化 设计 。 为 使 对 转 结 


的 


构 


攻 角 of? 8.068 6.008 4.118 2.726 1242 
1.3 喷 管 设计 方法 
喷 管 的 设计 通过 控制 过 流 面积 的 变化 来 获得 


均匀 稳定 的 流 场 ， 避 免 流 动 分 离 ， 减 小 损失 。 根 据 
收缩 比 1.77， 参 考 水 洞 收缩 段 设 计 经 验 
参照 文献 中 中 的 | 
积 变化 。 


定 喷 管 长 径 比 L/D=1。 
选用 五 次 方 曲 线 ， 将 其 转化 为 


H°, g 
| 线 数据 ， 
内 导 流 锥 


再 由 面积 变化 和 内 曲线 方程 计生 


曲线 为 抛物 线 ， 选 取 长 径 比 1， 得 到 其 
出 五 次 方 外 


WRITE. 
HRJ 


容易 实现 , MIKER 


K5, 从 而 简化 换 向 机 构 。 


考虑 到 前 后 转子 面向 来 流 的 位 置 不 同 ， 前 转子 更 侧 


男 外 ， 从 与 航行 体 的 匹配 性 能 > 


重 抗 空 化 性 能 设计 ， 后 转子 更 侧重 推进 性 能 设计 。 


看 ， 驱 动 系统 本 身 


程 ， 


fes 


EI 
S 


发 现 流 场 均 匀 稳 定 ， 无 流动 分 离 ， 说 明 喷 管 设 ? 


理 。 


并 对 分 界 点 进行 过 渡 光 滑 处 理 。 最 终 得 到 的 喷 
| 线 如 图 2 所 示 。 对 所 设计 的 喷 管 进行 数值 计算 ， 


人 
FT 


的 平稳 性 是 相当 重要 的 ， 为 消除 驱动 机 构 的 残余 


转 力矩 ， 前 后 转子 之 间 需 要 满足 扭矩 平衡 的 约束 条 


件 


避 型 设计 的 原则 在 于 辟 型 载荷 的 优化 。 考 虑 
转子 叶 项 部 分 空 化 性 能 变 差 的 问题 ， 利 用 辟 型 设 
冲 角 的 变化 实现 叶 项 辟 型 载荷 的 轻微 减少 ， 同 时 
到 削弱 叶 顶 泄露 流动 的 效果 。5 个 设计 刘 面 的 辟 


1H 
A 


到 
计 
达 


型 


均 采 用 NACA4 位 数 系列 ， 通 过 指定 经 型 骨 线 最 
点 的 相对 位 置 xp /1 从 叶 根 到 叶 顶 线性 增加 ， 即 载 
点 逐渐 后 移 的 方式 ， 进 一 步 实现 叶 顶 的 卸载 以 改 
其 空 化 性 能 。 表 3 为 最 终 确定 的 翼 型 设计 参数 。 


N 


表 3 EE 


Table 3 Design parameters of profile 


r 
Ted 
ET 
何 


E 
E: 
r1 


设计 截面 l1(Hub 2 3 4  S(Tip) 
o ifie 000 100 100 1.00 0.00 
h 相对 位 置 xy 0.50 0.525 0.55 0.575 0.60 
i 出 口 落后 角 5/ 10512 9.457 7.590 6.311 6.798 

JU ff al 8.617 5.902 3.328 1.984 0.863 
后 ”设计 冲 角 i/* X000 1.00 1.00 1.00 — 0.00 
$6 ”相对 位 置 xW 0.0 0.525 0.55 0.575 0.60 
CT 出 口 落后 角 5/” 13331 11332 9492 7731 7.165 
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Fig. 2 Nozzle curve 


1.4 数值 计算 方法 
本 研究 中 的 数值 计 旬 


M 
图 2 nie 


采用 商业 CFD 软件 


ANSYS CFX 16.0 进行 。 水 力 性 能 及 内 流 特 性 的 评 


估计 算 采 用 全 流 道 非 定常 计算 方法 ， 


估计 算 采 用 单 流 道 定常 两 相 流 计算 方法 。 


而 空 化 性 能 评 
采用 


道 网 格 数 约 52 万 ， 
人 泵 整体 单 流 道 总 网 格 数 约 为 87 73 


ANSYS TurboGrid 16.0 对 计算 模型 进行 网 格 划 分 ， 
) 关 性 分 析 ， 最 终 确定 的 前 转子 单 流 
后 转子 单 流 道 网 格 数 约 35 万 ， 


边界 条 件 设置 如 下 : 入 口 采 用 质量 流量 入 口 


(Bulk Mass Flow Rate), E H XA 
(Opening), "T Fr I $è Sx 3x He FE 
Wall), ^h5& 53K Hl sc I8] Jie fc E 


H FF AUN A A 
壁面 (Rotating 


(Counter Rotating 


Wall). 

空 化 计算 中 介质 选用 25 摄氏 度 的 水 和 水 蒸气 
共 两 相 ， 初 始 气相 体积 分 数 为 0， 初 始 液 相 体积 分 
数 为 1。 满 流 模型 采用 SST 模型 ， 空 化 模型 采用 
Rayleigh Plesset 模型 。 求 解 格式 为 高 阶 求解 模式 的 
差分 格式 ， 均 方 根 残 差 值 设 为 le-4。 


2 对 转 泵 喷 推 进 器 的 水 力 及 空 化 特性 


2.1 水 力 性 能 

对 转 和 泵 喷 推 进 器 的 非 定常 水 力 性 能 数值 计算 
结果 如 图 3 所 示 ， 设 计 工 况 下 泵 喷 扬 程 约 为 Ho— 
15.28 m, 泵 效率 w= 90%, 推力 T=1.48 kN。 由 图 3 
HJI, XPERT H-O 曲线 相对 较 陡峭 ,效率 曲线 
在 设计 工 况 附近 偏 平 稳 ， 这 与 对 转轴 流 泵 的 性 能 研 
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图 3 对 转 泵 喷 水 力 性 能 曲线 


Fig. 3 Performance curve of contra-rotating waterjet pump 


2.2 空 化 性 能 

为 得 到 研究 对 象 的 空 化 特性 曲线 ， 本 文通 过 改 
变 出 口 的 压力 值 进行 数值 计算 ， 得 到 一 系列 有 效 空 
化 余 量 NPSH, 和 总 扬程 五 的 对 应 关系 ,出 口 压 力 值 
从 150kPa 开始 一 直 降 到 50kPa， 对 应 NPSH, 值 从 
14.1m 降 至 7.7m。 每 次 计算 利用 上 一 次 的 计算 结果 
作为 初始 条 件 ， 以 提高 计算 效率 。 绘 制 出 的 对 转 泵 
喷 转 子 的 空 化 特性 曲线 如 图 4 所 示 。 


Him 


740 11 12 13 14 15 
NPSH m 
图 4 对 转 泵 喷 转 子 的 空 化 特性 曲线 


Fig. 4 Cavitation curve of contra-rotating waterjet pump 


取 扬 程 下 降 3% 时 的 NPSH, 作为 对 转 泵 的 空 
化 余 量 的 临界 值 。 由 空 化 特性 曲线 得 ， 临 界 空 化 余 
量 NPSH,=8.44m。. 传 统 的 空 化 比 转速 计算 式 如 式 (1) 
所 示 。 


5.62n JO 
C= DIr 
NPSH?* 


IP, n/min! 为 转子 转速 ，NPSH, /m 为 必需 
空 化 余 量 。 

按 式 (D) 来 计算 对 转 泵 的 空 化 比 转速 HX 
n=2000rpm， 用 NPSE 代替 NPSH,， 计 算得 空 化 比 
转速 C=1045。 

但 是 由 于 对 转 泵 采用 两 个 转子 对 旋 ， 用 传统 的 
空 化 比 转速 定义 式 不 能 体现 其 具体 的 空 化 性 能 。 日 
本 九州 大 学 流体 控制 研究 室 的 Furukawa APTA A 
推导 了 相同 设计 条 件 下 ， 对 转 式 轴 流 泵 和 传统 轴 流 
泵 所 满足 的 关系 式 如 式 (2) 所 示 。 


(1) 


Dy Npr 二 d (2) 


Disngs NE 


AP, 下 标 RR RS 分 别 代表 对 转 泵 和 单 转子 
Fo 

根据 式 (2)， 将 对 转 泵 换算 成 等 效 的 单 转子 泵 ， 
得 到 等 效 的 空 化 比 转速 C=1810。 

为 进一步 对 比 空 化 性 能 ， 我 们 选取 了 一 台 设 计 
条 件 相 同 的 单 转 子 泵 作 比 较 。 在 其 设计 工 况 下 ， 流 


t& O=0.212mz/s, {FE H = 13.15 m, RAZ y= 8995. 
计算 得 其 临界 空 化 余 量 NPSH,=12.6m, 空 化 比 转速 
C=1160。 由 此 可 见 对 转 泵 的 空 化 性 能 得 到 了 显著 的 


提高 。 


2.3 局 部 欧 拉 能 量 分 布 

取 叶 片 展 向 的 5 个 设计 截面 ， 研 究 局 部 欧 拉 能 
量 沿 着 轴线 方向 的 分 布 。 局 部 欧 拉 能 上 
究 流 面 从 转子 入 口 到 转子 出 口 的 能 量 增长 过 程 ， 其 
定义 式 如 式 (3) 所 示 。 


WP, m 代表 子午 方向 坐标 ，s 代表 展 向 坐标 ， 
U 为 圆周 速度 ， 乃 为 绝对 速度 的 圆周 分 量 ， 芭 为 绝 
对 速度 的 圆周 分 量 。 
对 转 和 泵 喷 和 同 设计 条 件 的 单 转子 泵 的 局 部 欧 


Pressure /kPa 


400 
0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 


拉 能 量 分 布 曲线 如 图 5 所 示 。 
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图 5 局 部 欧 拉 能 量 分 布 曲线 

Fig. 5 Local Euler head distribution 


Local Euler Head /m 


-0.02 0.00 


由 图 可 知 ， 对 转 泵 喷 的 两 个 转子 分 别 承担 了 一 
部 分 载荷 ， 而 单 转子 泵 则 只 靠 一 个 叶轮 承载 全 部 载 
fu. 另外， 前 转子 曲线 的 斜率 较 小 ， 说 明 叶 片 加 载 
缓慢 ， 后 转子 斜率 偏 大 ， 说 明 能 量 密度 高 ， 这 与 我 
们 等 速 差异 化 设计 的 理念 相 一 致 。 而 对 于 单 转子 泵 ， 
它 的 斜率 明显 更 加 陡峭 ， 说 明 叶 片 能 量 梯度 大 ， 更 
容易 发 生 空 化 和 二 次 流 现象 。 


2.4 叶片 载荷 分 布 

绘制 对 转 和 泵 喷 和 同 设计 条 件 的 单 转子 泵 的 叶 
片 载荷 分 布 曲线 如 图 6 所 示 。 由 图 可 知 ， 单 转子 泵 
上 叶 乒 的 最 低压 力 明显 低 于 对 转 泵 ， 而 且 载荷 分 布 
非常 不 均匀 ， 是 典型 的 前 部 加 载 。 而 对 转 泵 的 每 个 
叶片 上 的 最 低压 力 约 为 单 转 泵 的 一 半 ， 并 且 叶 片 载 
荷 基本 均匀 分 布 ， 为 中 后 部 加 载 。 因 此 对 转 泵 的 载 
荷 分 布 明 显 优 于 单 转子 泵 , 抗 空 化 性 能 得 到 了 提高 。 


Front rotor 
400 一 * 一 SPAN0.01 
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(a) Blade loading distribution of front rotor 
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(b) 后 转子 载荷 分 布 


(b) Blade loading distribution of rear rotor 


Streamwise (0-1) 
(c) 单 转子 泵 载荷 分 布 
(c) Blade loading distribution of RS type 


图 6 叶片 载荷 分 布 
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水 下 航行 体 性 能 参数 ， 通 过 相似 换算 设计 了 第 一 代 
实验 模型 对 转 泵 喷 及 其 实验 系统 。 简 述 了 对 转 泵 
转子 和 喷 管 的 设计 方法 ， 借 用 数值 计算 手段 ， 展 姑 
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Fig. 6 Blade loading distribution 


结论 


本 文 提 出 了 对 转 泵 喷 推 进 器 的 概念 ， 并 基于 茶 


pz 


mT 


对 转 泵 喷 推 进 器 的 水 力 性 能 和 空 化 特性 的 研究 ， 并 
与 相同 设计 条 件 的 一 台 单 转子 泵 进行 了 对 比 。 总 体 
而 言 结论 如 下 : 


1) 在 设计 工 况 下 ， 实 验 模型 对 转 泵 喷 扬 程 约 
X H= 15.28 m, FAK w= 90%， 推 力 7 
1.48kN X] £z BR] H-O 曲线 相对 较 陡 峭 ， 
效率 曲线 在 设计 工 况 附近 偏 平稳 。 

2) 实验 模型 对 转 泵 喷 的 临界 空 化 余 量 约 为 
NPS 太 =8.44m， 按 照 传统 空 化 比 转速 的 计 
算式 计算 其 空 化 比 转速 为 C=1045。 换算 成 
等 效 的 单 转子 泵 ， 得 到 等 效 空 化 比 转速 
C=1810， 大 于 同等 设计 条 件 的 单 转子 泵 。 
同时 ， 结 合 局 部 欧 拉 能 量 分 布 和 叶片 载荷 
分 布 可 知 ， 对 转 泵 喷 的 抗 空 化 性 能 得 到 了 


显著 的 提高 。 
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